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The propagation of fronts in the Fisher-Kolmogorov equation with spatially varying diffusion 
coefficients is studied. Using coordinate changes, WKB approximations, and multiple scales analysis, 
we provide an analytic framework that describes propagation of the front up to the minimum of 
the diffusion coefficient. We also present results showing the behavior of the front after it passes 
the minimum. In each case, we show that standard traveling coordinate frames do not properly 
describe front propagation. Lastly, we provide numerical simulations to support our analysis and to 
show, that around the minimum, the motion of the front is arrested on asymptotically significant 
timescales. 



I. INTRODUCTION 

Front propagation in reaction-diffusion equations 
(RDEs) is an important topic in the physical and biolog- 
ical sciences. These models appear throughout biology 
[TJ [2] , ecology [3] , cancer research dj , chemical kinetics 
[5], and geochemistry [6]. While classically concerned 
with homogeneous environments, interest in the behav- 
ior of fronts in heterogeneous environments has increased 
over the last several years. In particular, traveling waves, 
or propagating fronts, through certain classes of time [Jj 
and spatially varying environments [3J [E) have been well 
studied. Nonlinear, density-dependent diffusion, where 
the diffusion coefficient depends on u, has also been thor- 
oughly examined in p] [2] and [6] amongst others. The 
related problem of varying selectivity, or excitability, was 
examined in [8141 lj. For the reader interested in more 
mathematical issues concerning front propagation, we re- 
fer to P2J for a thorough review of this topic. Note that a 
complete listing of all references addressing these topics 
is beyond the scope of this paper, but the aforementioned 
references provide extensive bibliographies to the broader 
literature on front propagation in RDEs. 

We consider the case of a heterogeneous Fisher- 
Kolmogorov (FK), or Kolmogorov-Petrovskii-Piskunov 
(KPP), equation where heterogeneities are represented 
by spatially varying diffusion, i.e., 



u t = (a(x)u x ) x + f(u). 



(1) 



In the biological literature, spatially-heterogeneous diffu- 
sion coefficients a{x) are discussed repeatedly, but con- 
tinuously varying models are rarely investigated rigor- 
ously. For example, in ecological applications, several au- 
thors acknowledge the significance of heterogeneous dif- 
fusion [LMT7] and some have developed approaches for 
estimating a(x) across different habitats [T5JHH]. How- 
ever, many researchers dismiss spatially varying diffu- 
sion by attributing the diffusion variability to an evo- 
lutionary response in a sub-population (c/. P33 I2"U]). 
When researchers do allow diffusion to vary, it is typ- 
ically assumed to vary periodically in a square wave 
(c/. [3j HH I2T1J23] ) as would commonly be encountered 
in agricultural or urban contexts. In the spatially dis- 



crete, or lattice, context, researchers consider patches 
where diffusion is constant on each patch (c/. |24j and 
the references therein). 

A morhpogenesis phenomenon in which a(x) has been 
given serious consideration is in a general RDE model for 
the slug stage of Dictyostelium discoideum. In this case, 
diffusion is directly modulated by the spatial distribu- 
tion of a morphogen gradient concentration. The result- 
ing diffusion profile across a gap junction is a hyperbolic 
cosine, i.e., a(x) = cosh(a;). The RDE model with this 
form of a(x) has been investigated from an analytical 
(establishing a scale-invariant property for the generated 
wave pattern) [35] as well as computational [5S[ [57] per- 
spective. More recent work has considered the impact 
of letting a(x) = D + r\x 2 (for D constant and 77 small) 
in the context of the avascular growth phase of cancer 
[2"B] as well as on Turing bifurcations of standing wave 
solutions to RDEs [29] , 

None of these investigations, however, have generated 
analytical results for the impact on traveling waves of a 
continuously, strongly varying, diffusion coefficient. Thus 
the form of a(x) (motivated by the biologically justified 
choice of cosh(a:)) we choose to study is 



a{x) 



(2) 



From a modeling perspective, the diffusion coefficient 
represents an environment in which the necessary re- 
sources sustaining diffusion of the front dwindle and then 
grow as the front moves from left to right. Thus, it 
offers a good canonical model for a diffusion-mediated 
barrier. For example, one could use this mathemati- 
cal structure to model a geographic barrier such as a 
mountain. Further, this choice also allows us to study 
the effect of strongly varying diffusion, which is in con- 
trast to the case looked at in [TT]. This is because the 
choice in ([2 1 is more analytically tractable than the choice 
of hyperbo ic cosine while still maintaining the concave 
shape with a global minimum. In this work, we show, 
the strong variation creates several different asymptotic 
regimes through which the behavior of the front changes 
significantly. Lastly, it is in some sense an examination of 
the solution behavior on a subdomain of the periodically 
varying case studied in [3J. 
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Throughout the paper, we take as initial conditions for 
the step initial condition 



u(x, 0) = 



1, x<x c (0) 
0, x > x c (0) 



(3) 



where x c (Q) denotes the location of the step. It is well 
known in the traveling wave literature (|30|) that for this 
initial condition (and constant diffusion), the solution 
converges to a propagating front with speed determined 
by whether the nonlinearity induces pushed or pulled dy- 
namics. Initial conditions with slower rates of decay gen- 
erate a waveform with a faster wave speed. We label the 
location of the front via x c (t) with the convention that 
u(x c (t), t) = 1/2, so that the initial step is given by u = 1 
for x < x c (0) and u = for x > x c (0). We then look for 
fronts, i.e., solutions u(x,t) bounded between zero and 
one that propagate from left to right along the spatial 
coordinate. 

In the FK equation, the minimum of the diffusion co- 
efficient represents a turning point, i.e., the sign of a'(x) 
changes through the minimum. We develop, via multiple 
scales and WKB methods, an asymptotic description of 
the propagation of the front up to the point it crosses 
the turning point. Our analysis shows because of the 
choice of a strongly varying diffusion coefficient that the 
front does not propagate in a simple traveling coordi- 
nate system, i.e., u{x,t) cannot be written in the form 
u(x + ct). Instead, the front travels along curves of the 
form xe ct — xq, where Xq is some constant. This is to 
say that we can write u(x,t) = u(hi(x) + ct), and we 
can then think of the front as a traveling wave in a more 
generalized sense. 

As can be seen from the analysis in Section [TTJ the 
minimum of ([2| forces a fundamental shift in the behav- 
ior of the tail of u(x,t) as the tail crosses the turning 
point. This ultimately impacts the front by causing a 
shock like structure to form as it approaches the turning 
point. At the turning point, we show numerically that 
the propagation of the front is arrested on a timescale of 
0{\j^fe). This reflects the fact that at the turning point 
there is an asymptotically small amount of resources nec- 
essary for front propagation, and thus the front slows or 
is trapped. 

Beyond the turning point, we have preliminary ana- 
lytic results which explain the behavior of the front. We 
hypothesize that the front past the turning point is de- 
scribed by solutions to a stationary FK equation, in a 
traveling coordinate of the form In (a;) + ct. This again 
shows the need for a more general notion of traveling 
waves in the case of strongly varying diffusion. How- 
ever, more work is necessary to create a complete an- 
alytic treatment of the front beyond the turning point, 
especially for how the speed of the front is selected. This 
will be addressed in a later paper. We show via numer- 
ical simulation that our analytic treatment of the front 
up to the turning point is accurate. The simulations also 
show, as analytically predicted, that there is a fundamen- 
tal shift in the dynamics on either side of the turning 



point, and the numerics shows how fronts develop past 
the turning point. 

The structure of the paper is as follows. In Section [TT| 
we present our argument for how the front propagates up 
to the turning point. In Section III we present an argu- 



ment which shows how the FK equation smooths out step 
initial data. In Section [TV] we present a means of finding 
a traveling coordinate in which the FK equation is sta- 
tionary, and we then present a series of arguments which 
supports the idea that beyond the turning point, fronts 
evolve according to this stationary equation. Finally, in 
Section [Vj we present our numerical results. 



II. THE SOFT FRONT APPROXIMATION 

When x c (t) <C —1, we assume that the solution u has 
a soft front, i.e., a condition we define via the asymptotic 
relationship 

a(x)u xx < a'(x)u x . 

This asymptotic condition is denoted as the Soft Front 
Approximation (SFA). When the location of the front is 
such that |a; c (i)| ^ 1, we have that a(x) ~ x 2 . Thus, if 
for A > 0, u <~ l/x x in a region around the front, then 

a(x)u xx ~ A(A + l)x~ x , a'(x)u x ~ 2\x~ x , 

so that the soft front condition requires that < A <C 1. 
We see that this approximation is not valid for step or 
rapidly decaying initial conditions. This issue is ad- 
dressed in the next section, in which we demonstrate 
how the front softens from step initial conditions. As 
we will show, diffusion relaxes the rapidly decaying pro- 
file on short time scales and then makes the soft front 
approximation valid. 

Using the SFA, we now must solve the semi-linear hy- 
perbolic equation 

Ut ~ a'(x)u x + f(u). 

The method of characteristics then gives the system of 
differential equations 



dx 
~dt 



-a'(x), 



du 



/(«)• 



Using our choice of diffusion model a(x) = x 2 
choosing f(u) = u(l — u), we get the solution 



e, and 



u(xe 2( > t - to \t Q )e t - t ° 
U{Xl t] ~ 1 + uixe^-to) , i )(e*- t « - 1) ' (_1 



which follows the characteristics x(t) 



x e 



-2(t-t ) < j n 



the case that xq < 0, we see that information propagates 
towards the origin as desired for a traveling front. How- 
ever, if xq > 0, then all information again propagates to 
the origin, and thus the soft front model cannot describe 
a front propagating past the origin. Hence, we have a 
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mechanism that explains front propagation that follows 
a decreasing diffusion coefficient, i.e., a(x) is strictly de- 
creasing when x c {t) < 0. 

Therefore, given a snapshot of the profile u{x,t) at a 
time to where the SFA holds, we can describe the propa- 
gation of the front up to a neighborhood of the minimum 
of a(x). As can be seen from Q, the front does not 
propagate as a traveling wave, i.e., there is not a con- 
stant wave speed such that u(x, t) ~ u(x — ct). However, 
if we suppose that 

u(xe 2 ^'\t )^l, 

for xe 2 ^"* - 1 < x c (to), then we see that Q gives us that 
u(x,t) <~ 1. Likewise, if we suppose that 

u(xe 2{t - to) ,t a ) ~ 0, 

for a;e 2 '*~* ) > x c (t Q ), then we see that Q gives us that 
u{x,t) ~ 0. In this case then, taking logarithms, we get 
that the front travels along curves of the form 



ln(x(t))+2(t-t )=ln(x c (t )). 



(5) 



In this log-transformed spatial coordinate we can see the 
front propagates with speed c = 2. We refer to ([5| as 
a traveling wave coordinate (TWC), and we see that the 
TWC and the SFA allow us to compute a generalized 
notion of wave speed. This idea is studied further in Sec- 



tion IV where we extend our analysis beyond the turning 
point. 

We note that on a finite time scale we do not imagine 
that u(x, t) transitions globally to an algebraically de- 
caying profile. We are therefore arguing that the SFA 
holds in a region around x c (t) and this region must be 
matched, via intermediate layers, to the far field of u(x, t) 
which should have a much steeper decay profile. To begin 
to solve this problem, we look ahead of the front where 
< u <C 1. To study this regime, we linearize the FK 
equation around u = 0, i.e., we let u = ev, and then 
collect all terms in e. As is common, we assume that 
/(«) > 0, /(0) = /(l) = 0, /'(0) = 1, and /'(I) < 0, so 
that we get an equation for v(x,t) of the form 



v t = (a(x)v x ) x + v. 



We suppose v is given by the WKB ansatz, (c/. [9], [TT] . 
and [31] for examples of this approach) 

v{x,t) ^A{x,t)e^ x ^\ 

from which we get the leading order problem 

<t> t + a(x)4> 2 x + 1 = 0. (7) 

Using the method of characteristics and noting that ([7]) is 
a Hamilton- Jacobi equation (c/. 32 J, with Hamiltonian 
H(p, x) = 1 + a(x)p 2 , we get that 



dx 



±2Hyfa~(x), cf>(x, t) = (H 2 - l)t + cj)o{xo{x, t)), 



where H = y/H — 1 is a constant along the character- 
istics. The characteristics associated with the choice of 
diffusion coefficient ^ are given by 



x 



e e 



This expression is equivalent to 



e = \ xq- i/s„ + e 



and we therefore get that 



\Jx 2 + e j (x + \Jx 2 + e) 



•TO 



e e 



^2Ht 



^±2Ht 



+ Vx 2 + , 



Thus we see that for \x\ ^> y/e, the characteristics are 

to leading order given by xo ~ xe T2Ht . In general we 
see that information can propagate to or away from the 
origin. Therefore, ahead of the front, the decaying tail 
of u{x, t) propagates along characteristics of the form 
x ~ x e~ 2Ht for that portion of the tail between the 
front and the turning point. The portion of the tail be- 
yond the turning point follows characteristic curves of 
the form x ~ xoe 2Ht so that information travels away 
from the origin. Note, the value of H changes along 
different characteristics, and so we see that information 
propagates at different speeds along different character- 
istic curves. 

Through the turning point, or when |a;| <C we get 
the leading order behavior 



x Q - ^y/esmh(2Ht) + yx + ^-^j cosh(2Ht), 
so that for xq < we have 



x(t) ~ Ve ( -1+ Jl -2tanh(2£ft) + 2^|sech(2^)j 



and for xo > wc have 



(6) x(t) ~ yfe ( -1 + Jl + 2tanh(2ift) + 2-psech(2^t) ) 



Note, neither formula is useful on all time scales, but they 
are useful on times scales of 0(1). We see for xo < that 
the characteristics in the layer |x| <C y/e propagate away 
from the origin back to — yfe. Thus the characteristics in 
the outer region |x| 3> yfe collide in some intermediary 
layer with those in the inner region. The WKB analysis 
thus shows how the SFA eventually breaks down. 

Likewise, the WKB analysis shows that as the front 
approaches the origin something like a shock, or steep- 
ened front, must form. Determining this shock structure 
is beyond the scope of the current paper, though we are 
able to demonstrate it numerically in Section [Vj We can 
get some hint though of the dynamics of the front from 
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the following argument. Once the front has entered the 
inner region, i.e., |x c (i)| \/e, the FK equation becomes 

u t ~ eu xx + f(u). 

If the solution u follows the dynamics of a propagating 
front, then u is of the form 

fx — y/ect \ 
u[x,t) ~ u I — -p — I , 

so that we expect the front to be trapped in the inner re- 
gion around the turning point on timescales of 0(l/y/e). 
While our analysis at this time is not complete, the nu- 
merical results supports this hypothesis. 



III. SOFTENING SHARP INITIAL 
CONDITIONS 

In this section we explain how starting with step initial 
conditions pj), we can transition from the sharp to soft 
front regime. We note that in the life of the traveling 
wave, the softening discussed here happens before the 
SFA becomes valid. However, presenting it first would 
distract from the more important results in Section [IT] 
We present it now so that a complete description of the 
behavior of the wave before barrier transit is available for 
reference in Section llVl 

We assume throughout this section that x c (0) 3> 1. 
The asymptotic condition describing the sharp front is 
that 

a(x)u xx > a'(x)u x , 

which, if u ~ l/x x , implies that A > 1. We also choose 
a parameter L ^> x c 1 so that u x (±L,t) = 0. We 
choose Neumann boundary conditions to allow for ana- 
lytical tractability. We introduce the fast time T = - , so 
that, using the ansatz, 

u = uo(x, T, t) + eux(x, T,t) + ■■■ , 

we get the equations 

d T u Q = d x (a(x)d x u ) , 

d T ui = d x (a(x)d x u 1 ) + f(u ) - d t u Q . 

As for the leading order behavior uq, we need only solve 
the linear diffusion equation to find it. Using separation 
of variables in space and time, we write Uo(x, t,T) — 
4>{x;i)ip{T) which leads to the expansion for uq 



Sturm-Liouville theory (c/. [33]) ensures the functions 
cj) n form a complete, orthonormal set with respect to the 
norm 



|^>(a;)| 2 dx. 



Since we assume that a(x) > 0, then 



A„ / \(f> n (x)\ dx = 



(t>n{x)d x (a(x)d x <j) n ) dx 



a(x) \4> n ,x\ dx < 0, 



so that A„ < 0. Further, we see <j)Q is a constant corre- 
sponding to the eigenvalue A = 0. We set 4>o = \/l/2L 
so that 1 1 4>o 1 1 = 1. Likewise, using the step initial con- 
dition for u, we can find the initial condition for u Q (t) 
as 

*o(0) = " 



2L 

The remaining terms cr n (t) have the initial conditions 

cr„(0) = a( ^-d x (j) n {x c le). 

Moving to the second term u\{x, T, t), we see, using 
Duhamel's principle, that we can write Ui as 



(x,T,t)= / Y^ ln {t : s)4> n {x)e~\ x ^ T -^ds, 
J o n=0 



»i 



where 



ln(t,s) = j (-d t u + u (l-u ))4> n dx. (10) 

Since Ao = 0, we see that a possibility for a secularity, 
which means the asymptotic series becomes invalid on 
0(1) time scales (c/. [31]), arises from computing 



7o(£, s)ds, 



since if 7o(i, s) were independent of s then u 2 would have 
a term that growing linearly in T. Expanding the inte- 
grand in Equation ( 10 ) gives 



-d t Uo(x, t, s) + u (x, t, s)(l - u (x, t, s)) = 

0o(-«T(t) + O-o(t)(l-0OO-o(*))) 



u (x,T,t) = <T n (t)(/) n (x)e XnT ', 



(8) 



^((1 - 2cr ^o)o-n - o n )(\) n e 



|A„|s 



where <f> n and A„ solve the Sturm-Liouville problem 

d x (a(x)d x cj> n (x)) = X n (f> n , d x (j> n (±L) = 0. (9) 



n,j>0 
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Using the orthonormality of the functions <p n , then from 
above one has 

/ j (t,s)ds= T(-& (t) + a (t)(l - <fo<To(t))) 
Jo 

oo 2 

Thus, in order to remove the secularity, we enforce the 
condition 

(T = a Q (l- <f>o<ro), 
which has the solution 

1 

Mt) = 



+ ae~ 



where a = — </>o + 1/cto- 

As for the terms 7„, one has 



/ 7n (t, S )e-^ T - s ^d S = T((l - 2(To^o)a n - ( r n )e-l A "l T 
Jo 



E 

m,j>0 



|A m | — |A TO 



C e -(|Am| + |A,|)r _ e -|A„|T\ _ 

- A, \ / 



Note, | A n | — | A m | — | Xj | ^ due to the linear independence 
of the orthonormal eigenfunctions <p n (X). If we remove 
the terms of order Te~l A "l T as T — »■ oo, then we have 
that 

& n = (1 - 2cr ^o)o-„, 



or 



M0) £ 



CT " (T) (l + <7o(0)^o(e*-l)r 
Thus, on timescales T = C(l/e), or t = 0(1), one has 

u{x,T,t)=u (x,T,t)+O{e). 

From this, since for n > 1, one has by orthonormality 



l (x)dx = 0, 



it follows that the average of T, t) to leading order is 
given by cr (r)^ , or 



(u) - <7 (i)</> 



The quantity (u) is given by 



1 + 



L+x c 



(11) 



(u(;T,t)) 



2L 



This result shows that the average should increase ex- 
ponentially fast on time scales of 0(1) or less, which 
amounts to showing that the transition region of the so- 
lution u becomes smoother and "softer" rapidly. Having 
done this, the solution then enters the SFA regime. 



IV. FRONT BEYOND THE TURNING POINT: 
REDUCTION TO A STATIONARY EQUATION 

Since the SFA only leads to a model that describes 
propagation of a front up to the turning point, we must 
now find some other means of trying to describe propa- 
gation of the front past the turning point. To do this, 
by taking x ^> 1, so that a(x) ~ x 2 , we suppose that 
u = u(r/(x,t)), which means the FK equation becomes 



du 
dr\ 

By choosing 



d 2 u 

r it-iz = ( X7 ?*) 2 ^2 + ( x2 ^) x ~jz + /(")■ 



du 
d-q 



Vx = ±-, Vt = c, 
x 



we get the stationary FK equation 



d 2 u , .du „, . 

(-c±l)^+/(«) = 0, 



drj 2 

with coordinate 



r](x, t) = ± In \x\ + ct . 

The variable r\ is another instance of a (TWC), and so 
again we see that strongly varying diffusion requires a 
generalization of the definition of a traveling front. Lin- 
earizing around both the "+" and "— " cases from above, 
for the "+" case, we get roots to the characteristic equa- 
tion of the form 



A+(c) = 



(c-l)±((c-l) 2 -4) 1 /2 



and for the "— " case, we get 

A _ (c) = (c+l)±((c+l) 2 -4)V^ 

The front behaves like u ~ i]e Xri , when A is a double root 
om~ e Xr >, when A is a single root. In either case, we see 
that the decay rate of the front is not exponentially fast, 
but is instead only algebraically fast since 

e Xr > = |a;| ±A e Act . 

In order to ensure non-oscillatory decay, in the "+" case, 
one must take c < — 1, and in the "— " case, one must 
take c > 1. We define the characteristic curves of the 
TWC as the curves of constant r](x,t) which are given 

by 

± In \xo\ = ± In \x\ + ct, 

so that x(t) = xoe Tct . Since the characteristics are iden- 
tical in both of the relevant cases, we take the "— " case as 
our convention. We also see that, in contrast to the char- 
acteristic curves in the soft front case, that information 
is transported away from the origin. Thus, it appears the 
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TWC can provide a mechanism for propagation of fronts 
past the turning point. 

However, the algebraic decay is unexpected, and at 
first glance would seem to imply the TWC does not pro- 
vide relevant information about the propagation of true 
fronts, i.e., it is not clear that 



lim u(x,t) = u(r/(x,t)), 



(12) 



where u(x, t) is some solution to the FK equation with 
arbitrary initial condition and u{rf) denotes a solution to 
the stationary FK equation. It is non-trivial, even for the 
constant diffusion case, to show that solutions of the FK 
equation for some class of initial data satisfy (12). The 



literature on this issue is large, so we only refer the reader 
to the foundational papers [30l I35H37] as an introduction 
to this issue. Lacking any rigorous proof of (12), we 



provide a formal result to support the hypothesis that 
fronts evolve according to the TWC argument. 

Motivated then by classical approaches to the FK 
equation (cf. [I]), we again look at the linearized FK 
equation ([6]). Taking x 3> 1, so that a(x) ~ x 2 , separa- 
tion of variables can be used to solve ([6]). This gives the 
solution 



v(x, t) 



I \n(x)+ct 



where we take c > 3/4 so as to eliminate oscillatory so- 
lutions. Thus, we see algebraically decaying fronts are 
quite natural in this problem. This analysis supports the 
argument that solutions found via the TWC are in fact 
representative of true front behavior. Again, this is to 



say that ( 12 ) holds for some class of initial data. We also 



note the TWC fits into the framework of the WKB anal- 
ysis presented in Section [TT] since the characteristics past 
the origin found via the WKB analysis are the same as 
the TWC with r\ held constant. 




FIG. 1: Plot of Solution u(x, t) - L = 100, x c (0) = 

e = .1 



0.8 
0.6 
0.4 
0.2 



-35, 




-100 -50 



50 100 



FIG. 2: Contour Plot of Solution u(x,t) - L = 100, 
x c (0) = -35, e = .1 
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NUMERICAL RESULTS 



FIG. 3: Solution Profiles at Various Times - L = 100, 
af c (0) = -35, e = .1 



We have developed an explanation for how a strongly 
varying diffusion coefficient affects the propagation of a 
front up to the turning point of a(x). Using the TWC, 
we have some analytical results explaining front propa- 
gation past the turning point. In this section, we corrob- 
orate our results with numerical simulations. To simu- 
late solutions to the FK equation, we implemented the 
numerical scheme presented in |38j . an implicit/explicit 
method that is second order in space and first order in 
time. In these simulations we chose the nonlinearity 
f(u) = u(l—u), as well as Neumann boundary conditions 
u x (±L,t) = 0, and step initial condition as in (J3j) . 

Figures [T] and [2] depict the development of the solution 
in the (t,x) plane as a surface and as a contour plot, 
respectively. Figure [3] shows the front before, at, and 
after the turning point. In Figure |4j we compare the 
numerical propagation of the front and our asymptotic 



10r 
8 
6 
4 
2 



-50 



-Numerics 
-Asymptotics 



50 

x 



100 



FIG. 4: Comparison of SFA to Numerics - L = 100, 
af c (0) = -35, e = .1 



theory. The two curves represent the function x c>n (t) at 
which u num (x c<n (t),i) = 1/2, where u num (x,t) denotes 
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the numerical approximation to u[x,t). In the case of 
the asymptotic curve, we take to = 0, and use the same 
step initial condition as used in the numerics in @. The 
agreement between the SFA and numerics is convincing, 
and thus we have confirmation that the SFA is valid up 
to the turning point. As can be seen the SFA breaks 
down at the turning point. From Figures [T] [2j and[3j we 
see that the front steepens substantially. 

Once the front reaches the turning point, we see in 
Figures 111 [2] and [3] that the front becomes trapped. In 
Section [Hi we have predicted the timescale of trapping 
to be 0(l/y/e). In Figure [H] we plot the time duration 



^7r 
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21 
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Numerics 
-Least Squares: -0.131+1 .51 (1/e) 1/ 



10 15 
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FIG. 5: Trapping Times in Layer Around Minimum - 
L = 100, x c (0) = -35 

that the simulated front spends near x = as a function 
of 1/e, where we define near as \x\ < 0.4. Note the value 
0.4 is the width of the spatial mesh step in the numerics, 
and thus the smallest scale on which phenomena can be 
distinguished. Figure [5] shows that a least squares fit of 
the times to a curve growing like 1/y/e is accurate. Note 
that a decrease the spatial mesh step size will reduce 
the error between observed time and the least squares fit 
curve (results not presented here). 

On the other side of the turning point, as seen in Figure 
[2j the solution propagates in a manner consistent with 
the log-transformed TWC, i.e., information is propagat- 
ing along curves of the form r\ = — In |a;|+ci, or x = Xoe ct . 
However, it is currently an open question as to how c is se- 
lected (as opposed to the constant coefficient case where 
the regularity of the initial condition chooses c) and this 
is a direction of future work. 



VI. CONCLUSIONS AND FUTURE WORK 

We have shown in the paper that, after diffusion 
smooths the initial conditions such that SFA is valid, 
the choice of a strongly varying diffusion coefficient with 
a global minimum implies 

• The definition of a traveling front must be general- 
ized via the Traveling Wave Coordinate as defined 
m Section HU 

• Using WKB analysis, we show the minimum, or 



turning point, of a(x) causes the formation of shock- 
like behavior and leads to a trapping of the front 
on asymptotically long time scales. 

• The behavior of the front on either side of the turn- 
ing point is fundamentally different, and on either 
side, the TWC is necessary to describe dynamics. 

• After the front is past the turning point, the TWC 
allows us to transform into a stationary FK equa- 
tion, thus simplifying all subsequent analysis. 

To the best of our knowledge, the case of strongly vary- 
ing diffusion has not been studied in the literature, and 
this article is the first investigation of traveling waves 
with continuously and strongly varying spatial diffusion. 
We think these results will prove especially useful in mod- 
eling real world data for two reasons. 

First, quadratic spatial diffusion could be used to 
model barriers such as mountains and our results will give 
estimates of a time to barrier transit. For example, the 
invasion of the midwestern US by gypsy moths has been 
extensively studied (c/. [39] and references therein). It is 
also well known that the moths' spatial diffusion rate is 
tightly correlated with the local habitat [40] . By estimat- 
ing the population growth rates and fitting a quadratic 
curve to the local diffusion rates from capture-mark- 
recapture experiments (like those described in [H]), the 
method developed in this paper could be used to esti- 
mate the waiting time until the moths reach a certain 
habitat. In a non-ecological context such as the mor- 
phogen gradients discussed in these barriers could 
be nutrient deficient regions in gap junctions. Here, our 
results could be used to estimate the time till the con- 
centration of morphogen or metabolite at some boundary 
reaches a critical threshold. 

Second, with our work, we now have a non-trivial spa- 
tially varying FK equation that (via a traveling wave 
ansatz) could be reduced to a much simpler stationary 
problem. This will allow for a greater degree of control 
and flexibility in curve fitting routines as it will only re- 
quire a snapshot at a single point in time. With the 
gyspsy moth invasion example, one could envision es- 
timating the wave speed from a single snapshot. Fur- 
thermore, the authors currently have a manuscript in 
progress applying these tools to the gypsy moth inva- 
sion in an effort to estimate the time to barrier transit 
and invasion wave speed. 

Naturally, there are also many future directions for the 
analytical aspect of this work. It would be of particular 
interest to develop asymptotic matching techniques that 
would allow for connecting the various regimes of the 
front, i.e., the sharp to soft front transition, and then 
the trapping layer around the minimum of a(x). Like- 
wise, it would also be of interest to develop higher order 
expansions for the speed c in the SFA regime. 
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